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Abstract 

Background: Tlie widespread use of genome sequencing provided evidences for tlie liigli degree of 
conservation in innate immunity signalling pathways across animal phyla. However, the functioning and 
evolutionary history of immune-related genes remains unknown for most invertebrate species. A striking 
observation coming from the analysis of the pea aphid Acyrthosiphon pisum genome is the absence of 
important conserved genes known to be involved in the antimicrobial responses of other insects. This 
reduction in antibacterial immune defences is thought to be related to their long-term association with 
beneficial symbiotic bacteria and to facilitate symbiont maintenance. An additional possibility to avoid 
elimination of mutualistic symbionts is a fine-tuning of the host immune response. To explore this 
hypothesis we investigated the existence and potential involvement of immune regulators in aphid agonistic and 
antagonistic interactions. 

Results: In contrast to the limited antibacterial arsenal, we showed that the pea aphid Acyrthosiphon pisum 
expresses 5 members of Macrophage Migration Inhibitory Factors (ApMIF), known to be key regulators of the 
innate immune response. In silico searches for MIF members in insect genomes followed by phylogenetic 
reconstruction suggest that evolution of MIF genes in hemipteran species has been shaped both by differential 
losses and serial duplications, raising the question of the functional importance of these genes in aphid immune 
responses. Expression analyses of ApMIFs revealed reduced expression levels in the presence, or during the 
establishment of secondary symbionts. By contrast, /\pMFs expression levels significantly increased upon challenge 
with a parasitoid or a Gram-negative bacteria. This increased expression in the presence of a pathogen/parasitoid 
was reduced or missing, in the presence of facultative symbiotic bacteria. 

Conclusions: This work provides evidence that while aphid's antibacterial arsenal is reduced, other immune 
genes widely absent from insect genomes are present, diversified and differentially regulated during 
antagonistic or agonistic interactions. 
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Host-symbiont interactions, Acyrthosiphon pisum 



(bmc 

Genomics 



* Correspondence: Christine.coustau@sophia.inra.fr 

^Sophia Agrobiotech Institute, INRA-CNRS-UNS, UMR 7254, 400 Route des 

Chappes, 06 903 Sophia Antipolis, France 

Full list of author information is available at the end of the article 

O© 2014 Dubreuil et a!.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BIoIVIGCI C6ntr8l commons Attribution License (http://creativecommons.Org/licenses/by/4.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain 
Dedication waiver (http://creativecommons.0rg/publicdomain/zero/l.O/) applies to the data made available in this article, 
unless otherwise stated. 



Dubreuil et al. BMC Genomics 2014, 15:762 
http://www.bionnedcentral.conn/1471 -21 64/1 5/762 



Page 2 of 1 2 



Background 

Invertebrate immune systems have been shown to include 
highly diversified recognition systems, complex regulatory 
processes, and specific effectors [1-5]. However, the func- 
tioning and evolutionary history of immune -related genes 
remains unknown for most invertebrate species. 

For example, the recent sequencing of the first aphid 
genome {Acyrthosiphon pisum) provided new insights into 
insect immunity but also raised novel evolutionarily- and 
functionally- intriguing questions [6]. Indeed, genes par- 
ticipating in the immune responses of other insects are 
missing in aphids, such as the microbial recognition genes 
PGRPs (peptidoglycan receptor proteins) involved in the 
recognition of both Gram-negative and -positive bacteria 
[7,8]. Even more intriguing is the absence of most of the 
components of the immunodeficiency (IMD) signaling 
pathway, which is critical for fighting Gram-negative 
bacteria, and the absence of most of the antimicrobial 
peptides conserved in other insects [7,8]. These results, 
further supported by expression studies [8], showed 
that the aphid immune system greatly differs from that 
of other well studied insects and it was hypothesized 
that aphids compensate for a deficient immune system 
by symbiont-mediated host protection and an extraordin- 
ary reproduction rate [8]. Aphids are plant sap-feeding in- 
sects associated with the obligatory nutrient-providing 
symbiont Buchnera aphidicola [9,10]. Aphids also interact 
with several species of facultative bacterial endosymbionts 
[11] that are found free in the hemolymph [11,12]. Interest- 
ingly, secondary symbionts have been shown to influence 
important fitness-related traits such as body pigmentation 
[13], offspring production [14] and resistance to parasitoids 
[15-17] or pathogens [18]. It has been hypothesized that 
these intimate aphid-symbiont interactions may have led to 
the loss of specific antibacterial compounds or pathways to 
accommodate long-term host-symbiont coevolution and 
symbiont maintenance [19]. 

An additional possibility to avoid elimination of mutual- 
istic symbionts is a fine-tuning of the host immune re- 
sponse [20,21]. Aphid immune system may be particularly 
well regulated to prevent or limit damage to their bacterial 
mutualists under pathogen infections. As a first approach 



to explore this hypothesis, we searched A, pisum genome 
for immune regulators, and we noticed the existence of 
five genes coding for Macrophage Migration Inhibitory 
factors (MIFs) within the A. pisum genome. In vertebrates, 
MIFs are important pro-inflammatory cytokines acting on 
key ceUular processes of the immune response such as ceU 
proliferation and apoptosis [22,23]. MIFs were identified 
in a variety of species, including protozoan, nematode, 
mollusk and crustacean species [5,24-29]. It was shown 
that a mollusk MIF (BgMIFl) not only presented the 
expected activities on cell proliferation and apoptosis 
but played a major role in the response against parasitic 
infection [5]. This study provided functional evidence of 
the conservation of major immune -related functions of 
MIFs in an invertebrate and highlighted the importance 
of these immune regulators in invertebrate immunity [5]. 
Interestingly, all species investigated so far presented one 
or two MIF gene copies. The existence of additional MIF 
copies in A, pisum therefore requires dedicated evolution- 
ary and functional analyses to better understand the com- 
plex immune interactions of aphids with their symbionts 
and pathogens. 

Here we characterized A, pisum MIF members and ex- 
plored their existence and evolutionary history in various 
insect phyla. To gain insights into their functional role, 
we analyzed the expression of the several MIF copies dur- 
ing aphid-symbiont-pathogen interactions. 

Results 

Acyrthosiphon pisum MIF (ApMIF) family members 

The search for predicted MIF gene sequences in the anno- 
tated genome of Acyrthosiphon pisum [6] (Acyr_2.0 as- 
sembly) revealed 5 MIF members referred to as ApMIF 1 
to ApMIF5 (Table 1). Using BLAST searches [30], we 
found one or several full-length ESTs for each ApMIF 
member supporting their existence, structural annota- 
tion and expression. Complementary searches evi- 
denced an additional hit, corresponding to a 6th MIF 
gene that was actually annotated in a previous version 
of the genome assembly (Acyr_1.0) (Table 1). BLAST 
searches using the originally predicted ApMIF6 cDNA se- 
quence against the 214920 redundant ESTs resulted in a 



Table 1 Accession numbers and summary information on ApMIF genes and their products 





Aphidbase 
Genes 


NCBI 
Genes 


NCBI 
RefSeq mRNA 


NCBI 
RefSeq proteins 


PF01187 
e-value 


Exp. 
validation 


ApMIFl 


ACYPI002465 


LOCI 001 61 225 


NM_001 162635.1 


NP_001 156107.1 


1 .53e-38 


Yes 


ApMIF2 


ACYPI006088 


LOCI 001 651 24 


XM_001 946905.2 


XP_00 1946940.1 


5.85e-24 


Yes 


ApMIFS 


ACYPI000036 


LOCI 001 44890 


NM_001 126157.2 


NP_001 11 9629.1 


4.61 e-25 


Yes 


ApMIF4 


ACYPI003547 


LOCI 001 62394 


NM_001 162060.1 


NP_001 155532.1 


2.51e-24 


Yes 


ApMIFS 


ACYPI005907 


LOCI 001 64929 


XM_00 1948047.2 


XP_00 1948082.1 


6.43e-15 


Yes 


ApMIF6 


ACYPI002954 


LOCI 001 61 756 


Withdrawn 


Withdrawn 




No 



Experimental validation (Exp. validation) consisted in the amplification and resequencing of cDNAs. 
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partial and fragmented alignment with only 4 similar ESTs 
that did not support the existence of a full length tran- 
script from a 6th MIF gene. In order to further 
characterize ApMIF members, RT-PCR amplification of 
the complete coding sequences of the predicted ApMIF 
genes was performed from RNA of two genetically dif- 
ferent lineages of A. pisum (LLOl and YR2). Amplicons 
were obtained and resequenced for ApMIFl to -5. Se- 
quence results confirmed the expression of full length 
transcripts of ApMIFl to ApMIFS. Amplification of a 
full length cDNA from ApMIF6 (based on originally 
predicted sequence) was however not obtained. 

The five predicted ApMIF proteins showed highly sig- 
nificant hits with Macrophage Migration Inhibitory factor 
domain (PF01187/IPR001398) (Table 1). Interestingly, se- 
quences of the 5 ApMIF predicted proteins differ consid- 
erably among each other (Additional file 1: Figure SI) and 
present 33% to 55% similarity in their sequences. MIF pro- 
teins from mammals and snails have been shown to 
catalyze the ketoenol isomerization of small aromatic sub- 
strates such as hydroxyphenylpyruvate and L-dopachrome 
methyl ester [5,31]. Post-translational cleavage of the initi- 
ating methionine exposes an N-terminal catalytic proline 
(Pro2) that is essential for MIF tautomerase activity [32]. 
Three of the five predicted ApMIF proteins present this 
Proline2 residue as well as other conserved residues such 
as Lys33, Ile65 and Tyr93 and Vall07, participating to 
MIF active sites [33] (see Additional file 1: Figure SI for 
alignment of ApMIF sequences). Another conserved fea- 
ture of MIF proteins is to be secreted via non-classical 
pathways, involving, for example an ATP binding cassette 
transporter [34]. As expected, no signal peptides were de- 
tected in ApMIF sequences (SignalP 4.0; [35]) but the 



SecretomeP server [36] predicted ApMIF 1 and ApMIF 2 
as secreted proteins (secP score 0.934 and 0.78 respectively). 

MIF family members in insects 

In order to further characterize this multigenic family, 
we first searched for MIF genes in sequenced insect ge- 
nomes. Intriguingly, with the exception of two MIF genes 
in Tribolium castaneum and in Bombyx morU no MIF 
genes were identified from the other insect genomes avail- 
able at the time of the study (Table 2). We then performed 
an extensive search for MIF transcripts from 27 hemip- 
teran species (Table 3) with available ESTs in public data- 
bases. Interestingly, MIF transcripts were detected, in 
particular, in phloem-sap feeding species belonging to the 
Delphacidae, Aphididae, Pemphigidae, Pseudococcidae 
and Psyllidae families (Table 3). MIF genes were identified 
in all Aphididae with available ESTs with up to 5 mem- 
bers. Note that except for A. pisum the number of MIF 
members is probably underestimated due to the absence 
of annotated genomic data and to the modest size of some 
of the transcriptomic databases. Aphid genomic and tran- 
scriptomic data are currently growing rapidly and will 
allow a complete identification of aphid MIF members in 
a close future. 

Phylogenic reconstruction of MIFs from hemipterans 

Phylogenic reconstruction of complete predicted MIF 
sequences from hemipterans was performed using the 
crustacean Daphnia pulex as an outgroup (Figure 1). 
With the exception of A, pisum ApMIF5, all predicted 
MIF proteins from Delphacidae and Aphididae clustered 
into five clades including one clade specific of Delphaci- 
dae and 4 clades specific of Aphididae (Figure 1). Each 



Table 2 Survey of the number of MIF genes In Insect genomes 

Species Nb MIF 

Diptera Drosophila species (r5.46, ftp://ftp.flybase.org) 0 

Anopheles gambiae (r3.6, ftp://ftp.flybase.org) 0 

Aedes oegypti (rl.3, http://aaegypti.vectorbase.org/) 0 

Lepidoptero Bombyx mori (http://www.silkdb.org/) 2 

Coleoptera Tribolium castaneum (v2.0, http://www.hgsc.bcm.tmc.edu/) 2 

Hymenoptera Acromyrmex echinatior (v3.8, http://www.antgenomes.org/) 0 

Apis melifera (r4.5, http://hymenopteragenome.org/beebase/) 0 

Apis florea (vl.O, http://www.hgsc.bcm.tmc.edu/) 0 

Bombus impatiens (v2.0, http://hymenopteragenome.org/beebase/?q=gbrowse_bimp) 0 

Bombus terrestris (vll, http://www.hgscbcm.tmcedu/) 0 

Nasonia vitripennis (vO.5, http://www.hgscbcm.tmcedu/) 0 

Hemiptera Acyrthosiphon pisum (v2.l http://www.aphidbase.com) 5 

Anoplura Pediculus humanus (ul.2, https://www.vectorbase.org/) 0 

Cladocera Daphnia pulex (vl.O, http://genome.jgi.doe.gov/pages/blast.jsf?db=Dappul) 3 



Between parentheses are indicated the assembly version numbers and the URL used for genomic database searches. Daphnia pulex is added as an outgroup. 
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Table 3 Survey of complete MIF family members from 27 available hemiptera species 



Order Suborder Family 


Species 


Num EST 


Nb MIF 


Hemiptera Auchenorryncha Cicadellidae 


Grophocepholo otropunctoto 


6481 


0 




Homalodisca vitripennis 


20030 


0 




Oncometopia nigricans 


9056 


0 


Heteroptera Delphacidae 


Niloporvoto lugens 


118020 


2 




Peregrinus moidis fPm) 


20678 


2 


Cimicidae 


Cimex lectularius 


7129 


0 


Miridae 


Adelpliocoris lineolotus 


2915 


0 


Reduviidae 


Dipetalogaster maxinnus 


2671 


0 




Rhodnius prolixus 


16105 


0 




Triotonno brosiliensis 


2115 


0 




Tfiotomo infestans 


2564 


0 




Triatoma motogrossensis 


2230 


0 




Triotomo rubido 


1850 


0 


Lygaeidae 


Oncopeltus fosciotus 


1115 


0 


Sternorrhyncha Aleyrodidae 


Bemisia tabaci 


11923 


0 


Aphididae 


Aphis gossypii (Ag) 


88851 


2 




Rhopolosipliunn moidis (Rm) 


17649 


4 




Rhopalosiphum padi (Rp) 


17892 


3 




Toxoptero citricido (Tc) 


4304 


2 




Acyrtliosiplion l<ondoi (Ak) 


23040 


2 




Acyrtliosiplion pisum (Ap) 


214834 


5 




Myzus oscolonicus (Ma) 


22942 


5 




Myzus persicae (Mp) 


27728 


3 




Sitobion avenae (Sa) 


3478 


1 


Pemphigidae 


Pemphigus spyrothecoe 


18587 


0 


Pseudococcidae 


Maconellicoccus hirsutus (Mh) 


7669 


2 


Psyllidae 


Diaphorina citri (Dc) 


19598 


1 


The total number of redundant ESTs (Num EST) available in public databases at the time of the study is shown for each species, as well as the number of 
complete cDNA unique sequences predicted to code MIF proteins (presence of PF01187 domain confirmed by pfam_scan). Between parentheses are shown the 
species identification letters as used for the phylogenetic reconstruction. 

ApMIF member clustered into one of the four clades, hemocytes, we used an antibody raised against peptides 
together with MIF members from other aphid species of one of the well-expressed ApMIFs (ApMIFl). As 
(Figure 1), showing that several duplications occurred shown in the Figure 2C, ApMIFl protein localizes within 
before the diversification of aphids but after the split of hemocyte granules. 

Aphididae from the other hemipterans. Further studies on ApMIF expression were carried out 

on the 5 ApMIFs. However, because expression levels of 
ApMIFs constitutive expression ApMIF4 and ApMIFS remained at the limit of detection 
Expression of the 5 ApMIF transcripts was analysed by by real time quantitative RT-PCR in each experiment, 
real time RT-PCR on 12-day-old aphids from YR2 and they are not presented in the following result sections. 
LLOl free of any secondary symbionts [14,37]. In both 

lineages, transcripts of ApMIFl, ApMIF2 and ApMIFS, Increased ApMIF expression after immune challenges 

respectively clustered in clade A, C and D, were highly To explore the expected involvement of ApMIF genes in 
detected in whole bodies while ApMIF4 and ApMIFS aphid immune response, we examined their expression 
showed very weak expression levels (Figure 2). All five levels after immune challenges. We first used Aphidius 
ApMIFs transcripts were detected in hemocytes, sup- ervi, which is the main parasitoid of A, pisum in the field 
porting their potential involvement in immune processes [38]. The two clones of A. pisum showed a significant in- 
(Figure 2B). To better assess the presence of ApMIFs in crease in ApMIFl, ApMIF2 and ApMIFS transcript levels 
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Figure 1 Maximum lil<elyiiood tree siiowing thie piiylogenetic relationsiiip between full-length) hemiptera MIF unisequences (listed in 
Table 3). Genbank accession numbers of the sequences used for the reconstruction are shown next to the species identification letters 
(2 letters as defined in Table 3). Clades that are specific of delphacidae or aphididae (A-D) are highlighted in light and dark blue, respectively. 
Values at nodes are bootstrap proportions. Bar: 0,2 substitution/site. 



after exposure to female parasitoid (Figure 3). An increased 
ApMIF expression is also observed in aphids after injection 
with the Gram-negative bacteria E, coli (Figure 3). 

Decreased ApMIF expression during interactions with 
mutualistic symbionts 

We first compared the expression levels of ApMIF tran- 
scripts in the three lineages derived from YR2 and each 
harbouring a specific secondary symbiont (YR2-Hd, YR2- 
Ri and YR2-Ss). Aphids carrying mutualistic symbionts 
showed lower expression levels of ApMIFl, ApMIF2 and 
ApMIFS (Figure 4). In order to further study changes in 
ApMIF expression in presence of secondary symbionts, 
we measured ApMIF expression levels in LLOl following 
the injection of each of the three secondary symbionts. 
Since it was reported that injection of hemolymph from 
infected to uninfected A. pisum can establish a stable in- 
fection of secondary facultative symbionts [39,40], each 
secondary symbiont was manually collected from YR2 
lineages and injected into LLOl adults. ApMIF expres- 
sion was measured 1 and 4 days after injection. ApMIF 1, 



ApMIF2 and ApMIFS expression levels rapidly decreased 
after infection as compared to controls injected with 
Schneider medium (Figure 5). To ascertain that experi- 
mental infections with the facultative symbionts were 
successful, individuals from the Fl progeny of injected 
adults were microscopically examined. Presence of bac- 
teria in the hemolymph confirmed the establishment of 
these vertically transmitted symbionts, and analyses of 
ApMIFs expression levels confirmed that their decreased 
expression persisted in the Fl progeny (Figure 5). 
Altogether, these results showed that the estabUshment 
and the presence of the three major facultative symbi- 
onts of A. pisum correlate with a decreased expression 
of the three well-expressed ApMIF genes. 

ApMIF expression during both agonistic and antagonistic 
interactions 

Because some secondary symbionts can impact important 
aphid fitness-related traits such as resistance to parasitoids 
[15-17] or pathogens [18], it may be disadvantageous 
for aphids to mount a complete immune response in 
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APMIFl APMIF2 APMIF3 APMIF4 APMIF5 




APMIFl APMIF2 APMIF3 APMIF4 APMIF5 



Figure 2 Expression ratios (normalized to EF1 expression) of the five ApMIF transcripts in whole bodies (A) and hemocytes (B) from 
the A. pisum genetic lineages LL01 (dark blue) and YR2 (light blue). Each bar represents the mean expression ratio ± SD obtained from three 
independent experiments. (C) Immunolocalization of ApMIFl in hemocytes. Merger of fluorescent micrographs showing two hemocytes with 
labeled actin (green fluorescence), nuclei (purple/bleu fluorescence) and ApMIFl protein (orange fluorescence). Scale bar represents 10 |jm. 



the presence of both these mutuaUsts and a pathogen 
or parasitoid. Similarly, if the negative regulation of 
aphid immune response results from a direct effect of 
mutualists on host immune system, one may expect 
that this regulation would be maintained after an im- 
mune challenge and would limit the increase in ApMIF 
expression. To test this hypothesis, ApMIF expression 
levels were analyzed after exposure to parasitoids in 
the YR2 lineages carrying each of the three facultative 
symbionts. No significant change in gene expression 
was observed for ApMIFl (Figure 6A). Increased ex- 
pressions were observed for ApMIF2 and ApMIFS 
after exposure to the parasitoid A. ervi (Figure 6B,C), 
the highest expression levels being observed in the 
YR2-Ss lineage. However, all expression levels remained 
far below those of the YR2 lineage without mutualistic 
symbionts (Figure 6). 

Discussion 

Widespread lineage- specific gene duplications have been 
observed in the A, pisum genome [6]. The number of du- 
plications in this aphid appears greater than that of any 
other sequenced insect with over 2000 gene families that 
have undergone gene duplication [6]. However, these 
duplications were not reported to particularly affect 
immune-related genes, and, on the contrary, important 
and conserved genes of the antibacterial pathways are 
missing from A. pisum genome [6]. 

A specific search for MIF genes in genomes from insects 
revealed that except for Tribolium castaneum and 
Bombyx mori, no other genome contained MIF members. 




r 




APMIFl 



APMIF2 APMIFS 




APMIFl 



APMIF2 



APMIFS 



Figure 3 Relative expression (normalized to EF1 expression and 
to expression in control aphids (black)) of the three well-expressed 
ApMIFs transcripts, three days after immune challenge by 
Aphidius ervi (dark blue) or E coli (light blue) in the A. pisum 
lineages LL01 (A) and YR2 (B). Each bar represents the mean 
expression ± SD obtained from three independent experiments. 
Significant differences as compared with control are indicated by 
stars (Anova on row data: ***p < 0.0002). 
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APMIFl APMIF2 APMIF3 



Figure 4 Comparison of ApMIFs expression in YR2 aphids associated with R. insecticola (YR2-Ri), S. symbiotica (YR2-Ss) or H. defensa 
(YR2-Hd). Expression has been normalized to EFl expression and to expression in control aphids (black). Each bar represents the mean 
expression ± SD obtained from three independent experiments. Significant differences as compared with control are indicated by stars (Anova on 
row data: p < 0.0002; **: p < 0.002). 



A 0'06 



o 

■■g 0,04 




Dl D4 Fl 




ApMIFl ApMIF2 ApMIF3 



Figure 5 Effect of secondary symbiont injection on ApMIFs 
expression in LL01 aphids injected with R. insecticola (A), S. 
symbiotica (B) or H. defensa (C). Expression was measured at one 
(Dl) or four (D4) days post injection, as well as in offspring from the 
first generation (Fl), in symbiont injected (blue) and buffer injected 
(black) aphids. Significant differences as compared with control are 
indicated by stars (Anova on row data: '''''': p < 0.0002; p < 0.002). 



Conversely, analysis of hemipteran transcriptomic data- 
bases provided evidence for the existence of a MIF mul- 
tigenic family in Aphididae species as well as for the 
existence of at least 2 MIF members in Delphacidae spe- 
cies. Phylogenic reconstruction of complete predicted 
MIF sequences from hemipterans showed that sequences 
clustered into five clades including one clade specific of 
Delphacidae and 4 clades specific of Aphididae. Each 
Aphididae clade contained closely related MIF members 
from several aphid species showing that duplications 
occurred before the diversification of aphids. Altogether, 
these results showed that losses occurred in most insect 
phyla while both differential losses and duplications oc- 
curred within the hemiptera phylum. This differential loss 
versus duplication, together with the previously observed 
absence of important immune -related genes in A. pisum 
genome, raised the question of the functional significance 
of MIF members in Aphididae. 

Detection of the ApMIF transcripts in aphid whole 
bodies and hemocytes is consistent with the expected 
expression of these cytokines. Mammalian MIFs are 
expressed in various cell types, including all immune- 
relevant tissues (e.g. lung and digestive epithelia, skin) 
and macrophages [22]. In mammals, MIF genes are con- 
stitutively expressed and the protein is stored in granules 
of immune-relevant cells and rapidly secreted in extra- 
cellular fluids upon challenge [41,42]. Here we showed 
that ApMIF proteins are also detected in immune cell 
granules, suggesting that they are stored and secreted 
upon challenge. In addition to the expression localization, 
the significant increase in ApMIF expression after im- 
mune challenges support the role of these genes in pro- 
moting and regulating major early processes of the innate 
immune response [22]. MIFs are pleiotropic proteins 
involved in numerous cellular processes and are key 
cytokines controlling the response to both parasitic 
and bacterial infections [5,22,43,44]. The over-expression 
of ApMIFs was observed in two distinct aphid clones and 
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Figure 6 Expression levels of ApMIFI (A), ApMIF2 (B) and ApMIFS (C) in the YR2 lineage without symbiont (YR2) or associated with 
R. Insecticola (YR2-Ri), S. symblotica (YR2-Ss) or H. defensa (YR2-Hd), at 3 days post-exposure to the parasitoid A. ervl (Ae). Significant 
differences in expression as compared witli control apliids (ctl) not exposed to A. ervi are indicated by stars (Anova on row data: p < 0.0002; 
p < 0.002). 



after exposure to both parasitoids and pathogenic bacteria, 
supporting a general involvement of ApMIFs in aphids 
immune response, regardless of the aphid genetic back- 
ground or of the pathogen. 

Interestingly, a lower expression of ApMIF genes was 
observed either constitutively or after immune challenge 
in aphid lineages infected by one of the three tested sec- 
ondary symbionts (R. insecticola, S. symbiotica and H, 
defensa), A recent study showed that the number of A. 
pisum hemocytes is significantly smaller in YR2 lineages 
carrying these secondary symbionts as compared with 
the YR2 lineage without secondary symbionts [12]. The 
decreased expression levels of MIF genes that we ob- 
served in these YR2 lineages may therefore result partly 
from their smaller hemocyte number and partly from a 
down regulation of expression. 



These results are consistent with previous observations 
on the association between the weevil Sitophilus and the 
Sitophilus primary endosymbiont (SPE) [45]. In contrast 
to aphids, weevil immune system presents several families 
of conserved antimicrobial peptides [45] and it was shown 
that the antimicrobial peptide coleoptericin-A (ColA) in- 
hibits bacterial cell division and is essential in the regulation 
of endosymbiont number and location [19]. In addition, 
although expression of antimicrobial peptides (AMP) was 
similar in symbiotic (carrying SPE) and aposymbiotic lar- 
vae, a lower AMP gene expression was observed after im- 
mune challenge with E. coli in symbiotic insects [45]. 
Here, we observed a negative effect of symbionts on expres- 
sion of another type of immune genes in aphids, further 
supporting the importance of immune gene expression in 
regulating both agonistic and antagonistic interactions. 
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Conclusions 

The present study shows that members of the Macrophage 
Migration Inhibitory Factor (MIF) present a complex evolu- 
tionary history characterized by differential losses and du- 
plications across insect phyla. 

The 5 MIF members of A. pisum are functional, 
expressed in circulating immune cells and differently reg- 
ulated during a pathogenic or mutualistic interaction. 

This work provides evidence that while aphid's antibac- 
terial arsenal is reduced, immune regulators widely absent 
from insect genomes are present, diversified and differen- 
tially regulated during antagonistic or agonistic interac- 
tions. This supports the hypothesis of a fine-tuning of the 
immune response possibly accommodating both symbiont 
maintenance and response to aggressors. 

Methods 

Database searches and sequence analyses 

Protein sequences encoded by MIF genes in A, pisum 
genome were identified using BLAST programs [30] at 
the Aphidbase (available at http://www.aphidbase.com/) 
and National Centre for Biotechnology Information 
(NCBI) servers, providing a 6.2X-coverage genome as- 
sembly (Acyr_2.0 and Acyr_1.0) of A. pisum [6]. Other 
public databases used included the Reference Sequence 
(RefSeq) database [46] for A, pisum mRNA and protein 
sequences and the 214,920 A, pisum ESTs available at 
the time of the study in NCBI-dbEST. 

An exhaustive search of Hemiptera MIF transcripts in 
the NCBIs databases (NR, dbEST) was performed using 
the BLAST program [30]. The predicted translation of 
each full-length or partial sequence was checked for the 
presence of the conserved MIF domain (PFAM refer- 
ence number PF01187 and/or InterPro reference num- 
ber IPR001398) using pfam_scan software (available at 
http://pfam.sanger.ac.uk/search). Quality of sequences 
and their annotations was validated by multiple align- 
ments using Muscle [47] . Only full-length proteins were 
then used, re-aligned with Muscle and cleaned using 
Gblocks (version 0.91b, available at http://www.phyl- 
ogeny.fr). A maximum likelihood phylogenetic tree was 
generated using PHYML sofware (version 3) [48] and 
the WAG model of amino acids subsititution matrice. 
Bootstrap values were obtained from 1000 replicate 
samples using the proportion of invariable sites and 
gamma distribution parameters defined by PHYML. 

Aphid lineages 

Two distinct clones of A, pisum were used in this study, 
LLOl [37] and YR2 [14], both carrying only the obligate 
symbiont Buchnera aphidicola. For YR2, we also used 3 
different lineages derived from the same genetic back- 
ground but each differing by their respective secondary 
symbionts, being Regiella insecticola (U-type or PAUSS), 



Hamiltonella defensa (T-type or PABS), or Serratia sym- 
biotica (R-type or PASS) [14]. These lineages are referred 
to as YR2-Ri, YR2-Hd and YR2-Ss, respectively. 

Aphids were reared on Vicia faba under long photo- 
period (16 hours light/8 hours dark) and temperature 
(19°C) conditions to maintain continuous parthenogenetic 
reproduction [49]. Parthenogenetic wingless adults (12- 
15 days old) were used for all experiments except for 
parasitism by Aphidius ervi that was carried on 3 days- 
old larvae. 



Hemocyte collection 

Hemolymph was collected in a total volume of 200 (iL 
of aphid Schneider medium ( Sigma- Aldrich), as previously 
described [12]. Hemolymph was centrifuged at 200 g for 
5 minutes to eliminate putative embryos and hemocyte- 
containing supernatant was either used for total RNA 
extraction or further processed for immunolabeling. 

Immunolabeling of ApMIFI in A.pisum hemocytes 

Hemocytes were fixed in paraformaldehyde 4% in phos- 
phate buffer 0.1 M pH7.2 for 10 minutes and washed 3 
times in PBS under agitation. Cells were then permeabilized 
with Triton XlOO at 0.1% in PBS for 10 minutes and 
washed 3 times with PBS under agitation. Hemocytes were 
incubated with a custom-made antibody raised against 
two ApMIFl peptides (shown in the Additional file 1: 
Figure SI) (diluted 1/1000 in PBSlX-0.05% tritonXlOO) 
for 1 hour at room temperature. Control hemocytes 
were incubated with rabbit pre-immune serum (diluted 
1/1000 in PBSlX-0.05% Triton XlOO). After three 
washes in PBS, hemocytes were incubated for 1 hour at 
room temperature in anti-rabbit Alexa Fluor 488 second- 
ary antibody (diluted 1/500 in PBSlX-0.05% Triton XlOO) 
and anti-phalloidin-X5-FluoProbe (diluted 1/200 in 
PBS-0.05% Triton XlOO) to detect actin cytoskeleton. 
After 3 washes in PBS, nucleuses were marked by incu- 
bation in DAPI (1 (ig/mL) for 10 minutes. Hemocytes 
were mounted with anti-fading agent (Roti®-Mount 
Fluorcare, Carl Roth) and observed with a microscope 
Axioplan2 (Carl Zeiss) equipped for epifluorescence 
microscopy and differential interference contrast optics. 
Images were collected with a digital Axiocam (Carl Zeiss). 

Bacterial challenges 

Adult LLOl and YR2 aphids were injected with ap- 
proximately 300 E. coli (DH5 a) in insect Schneider 
medium (Sigma- Aldrich) using a nanoinject II nanoin- 
jector (Drummond Scientific Company) under a dis- 
secting microscope. 

To study changes in MIF gene expression of LLOl aphids 
during secondary symbiont infections, 12-days-old aphids 
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from the lineages YR2-Ri, YR2-Hd or YR2-Ss were pricked 
in the abdomen to collect hemolymph in a total volume of 
200 [iL aphid Schneider medium ( Sigma- Aldrich). Pelleted 
secondary symbionts were resuspended in 50 [iL of insect 
Schneider medium and counted under a microscope, 
using a counting chamber. Approximately 300 facultative 
secondary symbionts were injected in each adult LLOl 
aphid using a nanoinject II nanoinjector (Drummond Sci- 
entific Company) under a dissecting microscope. Control 
aphids were injected with sterile Schneider medium. 
Injected aphids were either frozen at -80°C at one and 
4 days post-injection, or allowed to deposit nymphs for 
2 days. 12-days-old adults from the first generation (Fl) of 
injected aphids were collected and stored at -80°C until 
RNA extraction. 

Production of aphids parasitized by Aphidius ervi 

Three days old aphid nymphs from LLOl, YR2, YR2-Ri, 
YR2-Hd or YR2-Ss were individually exposed to a single 
Aphidius ervi female wasp for 24 hr. Aphids were col- 
lected 3 days after parasitoid exposure and stored at -80°C 
until RNA extraction. For each sample, a set of aphids 
was maintained on plants under normal conditions to 
observe mummification and validate the success of para- 
sitoid infection. 

Real-time quantitative reverse transcriptase-PCR 
(RT-qPCR) 

Total RNA was isolated from homogenized tissues (whole 
bodies or hemocytes) using the RNeasy micro plus kit 
(Qiagen) and quantified with a nanodrop (Agilent). First- 
strand cDNA was generated from 500 ng RNA using 
iScript cDNA synthesis kit, according to standard proce- 
dures (Biorad, California, USA). Gene specific primers 
were designed with Primer3 software [50]. Real-time quan- 
titative PCR was performed on a DNA Engine 2 (MJ Re- 
search, Minnesota, USA) with qPCR MasterMix Plus for 
SYBR green I (Eurogentec, Belgium) using one internal ref- 
erence transcript (elongation factor 1-a, NCBI RefSeq 
XM_001951252.2). Running parameters were 95°C for 
10 min, followed by 40 cycles of 95°C for 30 s, 60°C for 
30 s, 68°C for 30 s. All amplifications were performed in 
triplicate assays. Signal intensity was measured at the end 
of each elongation phase and results were analyzed using 
the Opticon 3.1 software provided by MJ Research. To as- 
sess the specificity of the PCR amplification, a melting 
curve analysis of the amplicon was performed at the end of 
each reaction and a single peak was always observed. 
Standard curves were established with four serial dilutions 
of first-strand cDNAs, ranging from 1/10 to 1/10000. 
ApMIFs expression levels were normalized to expres- 
sion of the internal control elongation factor 1-a (EFl) 
using the comparative CT method (Applied biosys- 
tems, USA). 



Statistical analyses 

All experiments were repeated three times independ- 
ently, and each sample corresponds to pools of 10-15 
individuals. Differences in relative ApMIFs gene expres- 
sion levels were tested for statistical significance by one- 
way ANOVA and the tukey-Kramer test on row data 
(Software Prism v.5.0, GraphPad). Expression of the tar- 
get genes are presented either as an expression ratio 
normalized to EFl expression or as a relative expression 
normalized to EFl expression and to expression in the 
control sample. 

Availability of supporting data 

The data set supporting the phylogenetic analysis is avail- 
able in the TreeBase repository (http://purl.org/phylo/tree 
base/phylows/study/TB2:S16212). 
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